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Abstract 

We present quantum graphs with remarkably regular spectral char- 
acteristics. We call them regular quantum graphs. Although regu- 
lar quantum graphs are strongly chaotic in the classical limit, their 
quantum spectra are explicitly solvable in terms of periodic orbits. 
We present analytical solutions for the spectrum of regular quantum 
graphs in the form of explicit and exact periodic orbit expansions for 
each individual energy level. 

PACS: 05.45.Mt,03.65.Sq 

1 Introduction 



Consider a point particle moving along a network of bonds and vertices. 
Schematically, the network is represented by a graph T (see Fig. 1 for an 
example), which consists of Nb bonds and Ny vertices. The vertices are 
denoted by V,- t \ a bond connecting vertices i and j is denoted by B^. The set 
of bonds and vertices of V defines its geometry. We define a set of the bond 
potentials, Uij(k,x), where x and k are correspondingly the coordinate and 
the momentum of the particle on the bond B^. The vertices of T may be 
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equipped with 5-sources, etc. The geometry of T does not uniquely define 
the dynamics of a particle on T. In fact, since for any given geometry the 
graph may be "dressed" with arbitrary bond and vertex potentials, there 
exist infinitely many "dynamical realizations" of T. We call the set of bond 
and vertex potentials the "dynamical dressing" of the graph. Previously 
[H, 0, |3], f|, H , mainly the "bare-bond" graphs were studied, where the particle 
moves freely on the bonds. 

In this paper we focus on cases that have no turning points on the bonds, 
i.e. the energy of the particle is larger than all of the bond potentials, E > 
Uij(x,k), x G Bij. A simple way to implement this condition is to require 
that the system is scaling || [7|, ||, ||, pTOfl . This implies Uij(x, k) = Ay(x) k 2 , 
where the functions Ay (a;) are bounded for all x. In this paper we consider 
only simple cases where the functions Ay(x) are x independent constants, 

Uij(x, k) = Xijk 2 . (1) 

This is very similar to moving on a free graph except for substituting the 
bond lengths with the action lengths 

&ij = PijLij, (2) 



where Ly is the length of the bond .By, and /3y = J 1 — Ay. The scaling 
assumption (|l|) is not an oversimplification of the problem. Plenty of room is 
left for very interesting phenomena. Moreover, scaling quantum systems of 
this kind are the analogues of certain electromagnetic ray-splitting systems 
which have already been investigated experimentally in the laboratory J7|, [5], 

For all but the most trivial graphs, i.e. linear or circular graphs with 
vanishing bond and vertex potentials, the classical motion on a graph, inde- 
pendently of any particular dressing, is fully chaotic with positive topological 



entropy [11]. This means that the number of possible periodic orbits traced 
by the particle increases exponentially with their lengths. If no dynamical 
turning points are present, the topological entropy is independent of the dy- 
namical dressing and depends only on the geometry of the graph. Since at any 
vertex different from a "dead-end" vertex the classical particle has to choose 
randomly between several possibilities (reflection, transmission, branching), 
the particle's dynamical evolution resembles a stochastic, Markovian process. 
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Given their classical chaoticity it is surprising that the density of states 
of quantum graphs can be obtained exactly in terms of periodic orbit ex- 
pansion series |3|, f|, Furthermore, quantum graphs are considerably 
"more integrable" than all the previously known exactly solvable quantum 
systems. For example we will show below that for a certain class of quantum 
graphs — we call them regular quantum graphs — there exists an explicit 
and exact periodic orbit expansion for every quantum energy level. In other 
words, although the classical limit of regular quantum graphs is chaotic, each 
individual level of their spectra can be obtained exactly and explicitly via an 
analytical formula containing an explicit sum over the periodic orbits of the 
graph. To the authors' knowledge this is the first time that the spectrum of 
a quantum chaotic system is obtained both exactly and explicitly. 

The formal definition of regular quantum graphs is based on the properties 
of the spectral equation [|, |5| 

det[l - S(k)} = 0, (3) 

where S(k) is the scattering matrix of the graph ||. The modulus of the 
complex function (^) is a trigonometric polynomial of the form 

cos(S k - 7r 7o ) - $(£;) = 0, (4) 

where 

= ^2aiCos(Sik - 7T7i). (5) 

i 

and 

So = tE / M x ) dx ( 6 ) 

is the total reduced action length of the graph V and the constant frequencies 
Si < S naturally emerge as combinations of the reduced classical actions (0) . 
Under the scaling assumption, the coefficients aj, 70 and 7$ are constants. 
We now define regular quantum graphs. They satisfy 

E I ail < 1- (7) 



a: 



The motivation for this definition is the following: it allows us to solve (|4]) 
formally for the momentum eigenvalues k n , 

7r r 1 f arccos($(A; n )), for n + n even , s 

k n = 7T [n + fi + 70 - 7T < rAn w c 1 j j (8) 

S S I vr — arccos(<P (£;„)), for n + fi odd, 



3 



where /z is a fixed integer chosen such that ki is the first non-negative solution 
of (||). Because of (0) the second term of (||) assumes only values between 
u and tt/S q — u, where < u = arccos(a)/S , < n/2S . Thus, for regular 
graphs, the points 

- 7T 

kn = ^-(n + 7), n = l,2,..., 7 = 7o + /i, (9) 

>->0 

are guaranteed not to be roots of ([|) and serve as separators between roots 
number n and n + 1. 

Obviously the function (|9|) reflects the average behavior of the levels of 
the momentum. It is simply the inverted average staircase, k^ = Nik)" 1 . 
Geometrically the points (H) are the intersection points between the staircase 
function, N(k) = ~ k n ), and the average staircase N(k) resulting in 

the crossing condition 

N(k n ) = N{k n ) = n. (10) 

The crossing condition flT0| ) is illustrated in (Fig. 3). In case if there are no 
source potentials defined at the graph vertices, it can be shown that 

7 = -JV(0) = i (11) 

The existence of the separating points @ implies that the roots @ are 
confined to the "root zones", or "root intervals" I n = [k n -i, k n ], n = 1, 2, . . .. 
If a < C < 1 holds (C constant), equation (§) implies the existence of finite- 
width root-free "forbidden zones" R n = (k n — u,k n + u) surrounding every 
separating point k n , where no roots of (|3|) can be found. The roots of (0) 
can only be found in the "allowed zones" Z n = \k n -i + u, k n — u], which are 
subsets of the root intervals /„ = \k n -i,k n ]. For C — > 1 the width of the 
forbidden regions shrinks, u — > 0, and the allowed zones occupy the whole 
interval, Z n — > I„. 



Since 5*o is the largest action in (f|) and (|5|), it can be shown [|TJ] that 
k n is the only root in Z n . Therefore there is exactly one root k n inside of 
Z n C /„, and this root is bounded away from the separating points k n -i and 
k n by a finite interval of length 2u. 

The existence of the separating points (|9|) and the root-free zones R n are 
the key for obtaining an explicit and exact periodic orbit expansion for every 
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root of (Q). The starting point for obtaining the explicit expressions is the 
exact periodic orbit expansion for the density of states, 



p(A0=£*(*-*i)- (12) 

3=1 



As shown in JT|, ||, f|, || it can be written explicitly as 

-I oo 

p{k) = p{k) + - Re E ^ E ^ e luS ° pk - (13) 
" p i/=i 

Here p(&) is the average density of states, ^ is the repetition index, and S®, 
A p are correspondingly the reduced action and the weight factor of the prime 
periodic orbit labeled by p. In the scaling case, S® and A p are fc-independent 
constants fll2"| . Multiplying the density of states by k and integrating from 
k n -i to k n yields the value of the root contained between these separating 
points, 

rk n rk„ °° 

/_ p{k)kdk= _ J2 6 (k-kj)kdk = k n . (14) 

Performing the same procedure using the series expansion representation (|i~3|) 
and the crossing condition (0), we obtain 

7T 2 ^-^ 1 ^ Ap n . 
^n = ^^--2^^E^ sin(-z/^ p ) sm(7rz/u; p n), (15) 

Do 7T p D p y=1 Z/ Z 

where u p = Sp/So, and the A p J s are assumed to be real (no vertex potentials). 

Since all of the quantities on the right-hand side of QT5| ) are known, this 
formula provides an explicit representation of the roots k n of the spectral 
equation (|3|) in terms of the geometric and dynamical characteristics of the 
graph. To our knowledge, this is the first time that the energy levels of a 
chaotic system are expressed explicitly in terms of a periodic orbit expansion. 
Previously, explicit formulae for individual energy levels were known only for 
integrable systems. In the context of periodic orbit theory, the energy levels of 
integrable systems are given by the Einstein-Brillouin-Keller (EBK) formula 
[|j~T|j . However, apart from a few exceptional cases Jl3|] EBK quantization is 
only of semiclassical accuracy. 
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The difference between ( |T3"D and ([15]) is profound. The density of states 



13) allows the computation of spectral points only indirectly as the singu- 



larities of (|T3"D. Formula flT5p, on the other hand, allows the computation of 
every quantum level individually, explicitly and exactly in terms of classical 
parameters. 

In order to demonstrate that the class of regular quantum graphs is not 
empty we present an explicit example: the one-dimensional scaled step poten- 
tial with Vo = XE. A sketch of this potential is shown in Fig. 3. Physically 
this potential is realized, e.g., by a rectangular microwave cavity partially 
loaded with a dielectric substance J7], [8|, [|. The scaling step potential is 
equivalent to the scaling three- vertex linear graph shown in Fig. 2 (b). It 
has two bonds L\ = b and L 2 = (3(1 — b); the single scaling constant (3 (see 
(|D) is given by (3 = y/l — X. The spectral equation is given by 

|det [1 - S(k)\ | = sin(LA;) - r sin [(^ - L 2 ) k] = 0, (16) 

where L = Li+L 2 , and r = (1 — (3)/ (1 + 0) is the reflection coefficient at the 
vertex V 2 between the two bonds. It defines the eigenvalues k n only implicitly 
and is usually solved by graphical or numerical methods. Application of (|T5|), 
however, solves (|T6|) explicitly in terms of periodic orbits such as the ones 
shown in (Fig. 2). In order to apply fllBD we need the coefficients A p . They 
are given by JTD|, [12] 

A p = (-1)*(»V(P)(1 - r 2 ) r ^/ 2 , (17) 

where r is the reflection coefficient at the middle vertex and a(p) and r(p) are 
correspondingly the number of the reflections and the transmissions through 
it. Since the reflection coefficient may be positive or negative depending 
on whether the particle scatters from the right or from the left, the factor 
(— l)xO) j s needed to keep track of how many times it appears with a minus 
sign, including the sign changes due to the wall (x — and x — 1) reflections. 

In order to illustrate the convergence of ( |I5|) we computed k±, k w and 
fcioo of the scaling step potential including periodic orbits of increasing binary 
length q. For the parameters of the potential we chose b = 0.3 (see Fig. 2) 
and A = 1/2. Figure 4 shows the relative error e$ = \k$ — fc^ xact | //c^ xact for 
n — 1, 10, 100 and q ranging from 1 to 150. We see that even for small q the 
relative error is very small, decreasing further for large q as a power-law in 
q. The power of convergence appears to be the same for all three k and is 
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close to —2. The convergence with q is an important result. It indicates that 
although (15) is only conditionally convergent it (i) converges to the correct 
result and (ii) is not just asymptotically convergent, but keeps converging 
when more and more periodic orbits are included. 

Additional examples of regular quantum graphs are provided by all linear 
and circular quantum graphs with at most two bonds per vertex, indepen- 
dently of the number of vertices. In other words, for any simply connected 
quantum graph and any dynamical dressing there always exists a set of scal- 
ing constants \j of finite measure such that the regularity condition (0) is 
fulfilled. Well-known particular cases of these simply connected quantum 
graphs are the "Manhattan potentials" , which are obvious generalizations of 
the simple step potential shown in Fig. 2 (a) to arbitrarily many steps inside 
of the well, and linear chain graphs with scaling 5 function potentials at the 
vertices. 

It should be emphasized that the "inverse staircase expansion" ([TJ) is not 
just a curious finding, valid for some simple ID systems such as quantum 
graphs. Similar explicit series may be obtained for more complicated higher 
dimensional systems when the following two key ingredients are available. 
The first ingredient is the exact series expansion of the density of states flT3|) , 
which has already been established for other classically chaotic systems such 
as, e.g., quantum billiards [JR]. The second ingredient is a (piercing) average 
staircase function N(k) or the inverted staircase function k n which intersects 
every stair of the staircase, N(k n ) = N{k n ) = n, n = 1, 2, .... The intersection 
points k n then serve as the separators for the possible root locations, and 
the procedure outlined in the text can be used to find the periodic orbit 
expansions for individual roots of the system at hand. In most cases, of 
course, it is highly nontrivial to obtain these two necessary ingredients. The 
quantum graphs themselves are an excellent illustration of this point. While 
the expansion ( |13|) is valid for all quantum graphs, it is the crossing condition 
( [TDD th & t is violated when the inequality ([?]) brakes down. The regular graphs 
are precisely those for which the line N(k) = Sq Ar /7r + N(0) satisfies (|T0|) and 
allows the application of the analytical procedure that resulted in the explicit 
formula ([L5D for the representation and computation of individual eigenvalues 
k n . 

Y.D. and R.B. gratefully acknowledge financial support by NSF grants 
No. PHY-9900730 and PHY-9984075; Y.D. and R.J by NSF grant No. PHY- 
9900746. 
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Figure Captions 



Fig. 1: A generic (quantum) graph with five vertices and six bonds. 

Fig. 2: (a) Simple step potential, a basic problem in one-dimensional quan- 
tum mechanics. Also shown are examples of Newtonian ("N") and non- 
Newtonian ( "NN" ) periodic orbits used in the periodic orbit expansion of its 
energy eigenvalues (see Text), (b) Three- vertex hydra graph corresponding 
to the step potential above. 

Fig. 3: The staircase function N(k) and the average staircase N(k). For 
the regular graphs the average staircase intersects every "stair" of the N(k) 
graph, with the separators (][]) showing as the intersection points N(k n ) = 
N(k n ). 

Fig. 4: Relative error e£> = \k$ - /4 exact) |//4 exact) of (H) (see text) by 
including periodic orbits up to length q. The three curves shown correspond 
to ki, k w and k wo as indicated in the figure. 
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